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Abstract. We present recent lattice results on QCD thermodynamics at non- 
vanishing baryon number density obtained from a 6th order Taylor expansion in 
the chemical potential. Results for bulk thermodynamic observables, in particular 
for fluctuations in the baryon number density, are found to be well described by a 
hadron resonance gas model at low temperature and an ideal quark gluon gas at 
high temperature. We also analyze the radius of convergence of the Taylor series 
and discuss the information it provides on the occurrence of a second order phase 
transition point in the QCD phase diagram. 



1. Introduction 

Two quite different aspects of the studies of QCD thermodynamics on the lattice 
gained most attention in the heavy ion community during recent years - studies of 
the QCD phase diagram at non-zero baryon chemical potential and the analysis of the 
influence of a thermal heat bath on basic properties of hadrons, eg. their masses and 
widths. Here we will concentrate on the former. For a recent discussion of in-medium 
properties of hadrons we refer to the proceedings of this years Quark Matter and 
Lattice |2] conferences. 

Lattice calculations at vanishing baryon chemical potential, /is, suggest that for 
physical values of the quark masses the transition to the high temperature phase 
of QCD is just a rapid crossover rather than a phase transition which would be 
signaled by singularities in bulk thermodynamic observables. On the other hand 
many QCD motivated model calculations suggest that for /is > and small values of 
the temperature the transition indeed is a first order phase transition. This suggests 
that somewhere in the interior of the phase diagram there exists a second order phase 
transition point (chiral critical point) as an end point of a line of first order transitions. 
For vanishing light quark masses this point would be a tri-critical point Zt which, with 
increasing value of the strange quark mass, is expected to move towards lower values 
of /is and eventually could reach the fiB — axis. 

The expected generic phase diagram for physical values of the quark masses is 
shown in Fig. ^ Its exploration for large temperatures and small values of the baryon 
chemical potential became accessible to lattice calculations through the application 
of different techniques such as Ferrenberg-Swendsen reweighting 4 , Taylor series 
expansions E] as well as through simulations with an imaginary chemical potential 
13 IHj . This led to a first analysis of the baryonic contribution to bulk thermodynamic 
observables like the pressure, baryon number density and various susceptibilities 
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Figure 1. Generic pliase diagram of QCD with pliysically realized quark mass 
values (left). The occurrence of a second order endpoint in the interior of the 
phase diagram is expected to be a consequence of the existence of a regime of 
first order transitions in 3-flavor QCD at vanishing fis and for small values of 
the up, down and strange quark masses. The chiral critical point then lies on a 
hyper-surface of second order transition points emerging from the fig = plane 
(right). 

pi 1101 rm I12j . Moreover, it provided first estimates for the location of the chiral 
critical point @] at which the transition to the high temperature phase of QCD turns 
from a rapid crossover into a first order phase transition; recent calculations suggest 
a critical value /x^'* ~ 400 MeV We will discuss these results in more detail 

in section 3. However, before turning to this it seems appropriate to briefly review 
results on the QCD transition at vanishing baryon chemical potential where we try to 
emphasize the role of the strange quark in QCD thermodynamics. We will do so in 
section 2. Section 4 is devoted to a comparison of lattice results on thermodynamics 
in the low temperature hadronic phase at fis > with a phenomenological approach 
in terms of a hadron resonance gas model. This has been discussed in more detail at 
this conference by K. Redlich 

2. The thermal QCD transition at vanishing baryon chemical potential 

For fj,B — the transition from the low temperature hadronic phase to the high 
temperature phase has been analyzed in many numerical calculations. In particular, 
it has been found that for the case of QCD with three degenerate light quarks the 
transition is first order. In fact, the use of improved discretization schemes for 
the fermionic part of the QCD Lagrangian has shifted this regime to rather small 
values of the quark masses which correspond to a pion mass even below its physically 
realized value jTHj. This strongly suggests that the transition is continuous for the 
physically realized spectrum of two light (m, (i)-quarks and a heavier s-quark. A second 
order phase transition can then only show up at non-vanishing chemical potential if 
the hyper-surface of 2"'^-order phase transition, in which also the line of 2"'^-order 
transitions found at /i^ = in the {ud, s)-plane lies, bends over the physical point. 
This situation is illustrated in Fig. fright). The left hand part of this figure shows 
the resulting phase diagram of QCD with the physically realized mass spectrum. 

The transition to the high temperature phase of QCD commonly is related to the 
restoration of chiral symmetry as well as deconfinement, i.e. the disappearance of an 
almost massless pion as a Goldstone particle and the sudden liberation of quark and 
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Figure 2. The light quark chiral condensate in QCD with 2 hght up, down and a 
heavier strange quark mass (open symbols) and in 3-flavor QCD with degenerate 
quark masses (full symbols) 1171 . The right hand part of the figure shows the 
pressure calculated in QCD with different number of fiavors as well as in a pure 
gauge theory 1181 . Note that (2+l)-flavor QCD here refers to QCD with two light 
quarks and a heavier (strange) quark with a mass proportional to the temperature, 



gluon degrees of freedom. This is, indeed, reflected in the temperature dependence 
of the chiral condensate jJT] and the the pressure ^Hl shown in Fig. |21 These figures 
also indicate that the presence of strange quarks with a mass rus ~ 0{T) has little 
influence on the thermodynamics in the low temperature phase. The temperature and 
quark mass dependence of the light quark chiral condensate in QCD with two light and 
a heavier strange quark is similar to that of QCD with 3 degenerate (light) quarks. At 
Tc, however, the strange quarks with a mass rris ^ do not yet contribute significantly 
to the pressure or energy density. Consequently the transition temperature in 2 and 
(2+l)-flavor QCD has been found to be quite similar. In the high temperature phase, 
on the other hand, bulk thermodynamic observables clearly reflect the presence of the 
additional strange quark degrees of freedom] |. 

This suggests that at least at = the strange quark has little influence on 
the dynamics of the transition. The calculations performed with improved actions 
| 17II18[[T^ indicate that for the physically realized quark mass spectrum the transition 
to the high temperature phase is signaled by a rapid but continuous change in bulk 
thermodynamic observables. The transition temperature of (2+l)-flavor QCD found 
in a recent analysis with the Naik action |17jf . 

Tc = (172 ± 11 ± 7) MeV , (1) 

is consistent with earlier findings for 2-flavor QCD based on simulations with the 
p4-action, = (173 ± 8 ± sys. err.) MeV [TS|. 

I In the calculation of the pressure in (2+l)-flavor QCD shown in Fig.|^the strange quark mass has 
been taken to be proportional to the temperature, i.e. ms/T =const. If one keeps instead nis fixed 
to its physical value the pressure will gradually approach the high temperature limit of 3-flavor QCD. 
This can be inferred, for instance, from the calculation performed in (20 . 

^ Here we have averaged over the two fits performed in I17J to extrapolate to the chiral limit and 
included systematic uncertainties arising from these different fits in the systematic error. 
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Figure 3. Temperature dependence of expansion coefficients for p/T"^ in 2-flavor 
QCD and for quark masses corresponding to a pseudo-scalar (pion) mass of about 
770 MeV. 



3. Taylor expansion of the pressure at non-zero baryon chemical potential 

As mentioned in tlie introduction different approaches have been foUowed to study the 
phase structure of QCD at non- vanishing baryon chemical potential and to determine 
basic thermodynamic observables. All these approaches are currently limited to the 
regime of large temperatures and fiq/T<l, where /i^ = ^b/^ is the quark chemical 
potential. We will concentrate here on a discussion of Taylor expansions of the 
partition function of 2-flavor QCD around /is = 0. 

At fixed temperature and small values of the chemical potential the pressure may 
be expanded in a Taylor series around /ig = 0, 

oo 

|r = ^i-2-E^"(r)(^)" , (2) 

where the expansion coefficients are given in terms of derivatives of In Z(V, T, /i^), 
1 9" In Z 

i.e. c„(T) = ~'q(^ — /T)"' series is even in (fj,q/T) which reflects the invariance 

of Z(V,T, fj,q) under exchange of particles and anti-particles.. The first coefficient, 
Co, simply gives the pressure at /x^ = shown in Fig. [21 right). The next non-zero 
coefficient, C2, is proportional to the quark number susceptibility at = j21| . 

Its rise with temperature has been attributed to changes in the interaction among 
quarks and anti-quarks in the vector channel |22| . Further Taylor expansion 
coefficients entering the calculation of the pressure have been analyzed now up to 
Oin^) [101123 El- The coefficients for n = 2, 4 and 6 are shown in Fig. [3 

In the high temperature, ideal gas limit, only the first three expansion coefficients 
are non-zero, co{T = oo) = 7nt7r^/60, C2(T = oo) = n[/2 and C4{T — oo) = nf/Air'^. 
As can be seen in Fig.[3|at T ~ 2To all coefficients reach their corresponding ideal gas 
values within ~ 20%. 

The next-to-leading order coefficient, C4, is strictly positive and develops a 
pronounced peak at the {fiq — 0) transition temperature, Tp, which in turn, also 
leads to a large peak in the quark number susceptibility once fiq/T>0.5 |1U| (see 
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Fig. 5(right)). This has been interpreted as an indication for an approach to the 
chiral critical point at which Xq is expected to diverge"*". 

A new feature shows up in the expansion coefhcients at 0{ii^). The coefficient 
cg is positive only below Tq and changes sign in its vicinity. This has important 
consequences for the analytic structure of the QCD partition function and gives first 
hints at the convergence properties of the Taylor series: If there exists a 2"'' order 
phase transition point in the QCD phase diagram this would be related to a Lee- Yang 
zero of the QCD partition function for some real and positive value of (fiq/T)'^. This 
critical value will be directly related to the radius of convergence, p{T), of the Taylor 
expansion, if the partition function has no further zeroes in the complex (/Xg/T)-plane 
closer to the origin. A sufficient condition for the radius of convergence to be due to a 
zero on the real (/i(j/T)-axis is that all expansion coefficients in the Taylor series are 
positive. This apparently is the case for all c„(T) analyzed so far at T<0.95Tq. 

Ratios of subsequent expansion coefficients provide an estimate for the radius of 
convergence of the Taylor series, 



These ratios are shown in Fig. 2[left). While the first ratio, C4/C2 is well determine 
in current lattice calculations the ratio cq/c^ still has large errors. It is, however, 
apparent that these ratios change drastically across Tq and seem to be approximately 
temperature independent in the low temperature phase. As we will discuss in the next 
section this is expected to be the case also for a non-interacting resonance gas which 
has an infinite radius of convergence. At present, however, the large errors on Pn{T) 
prohibit to draw a firm conclusion on the asymptotic behavior of p{T). 

Similar caution is requested when interpreting the large fluctuations in the quark 
number densities signaled by the rapid rise of the quark number susceptibility (see 
Fig. 5). As mentioned above we expect Xg to diverge at a 2"'' order transition point. 
The rapid increase of Xq is, however, partly due to the rapid rise of rig itself as well as 
due to the rapid increase in the pressure for T ^ Tq. A quantity reflecting the relative 
magnitude of fluctuations is given by the ratio. 



This quantity only depends on ratios of Taylor expansion coefficients, and should 
vanish at a 2"*^ order phase transition point. Its temperature dependence is shown 
in Fig. El^right). As can be seen the ratio rises rapidly across Tq and approaches the 
ideal gas value. It, however, does not show any sign for a drop in the vicinity of Tq 
that could be taken as evidence for the existence of a second order transition. 

4. Resonance gas versus lattice results 

Aside from the obviously rapid change of bulk thermodynamic observables at 
temperatures close to Tq it is a striking feature of the Taylor expansion that ratios 
of the expansion coefficients are, to a first approximation, temperature independent 

Fluctuations in the net quark number density are related to the isothermal compressibility, 
kt = —V~^ {dV/dp)j, = —Nq^ {dV/d^iq)rp ^ = Xq/'^qi which diverges at a 2"'* order transition 




(4) 




point. 
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Figure 4. Ratios of subsequent expansion coefficients of the pressure (left) and 
the dimensionless ratio Ap/T^Xq ft"' values of the chemical potential \iqlT 
(right). Horizontal lines show the expected results for a hadronic resonance gas 
and an ideal quark-antiquark gas below and above To, respectively. 



below Tq and distinctively different from the ratios above Tp. While for T > Tq 
the ratios rapidly approach the corresponding ideal gas values their low temperature 
behavior can be understood in terms of the properties of a hadronic resonance gas 
|25|. For T < Tq the baryonic contribution to the pressure arises entirely from 
hadrons, which have masses significantly larger than the relevant temperatures, i.e. 
Tribaryon ^ mp > SiTq. Their contribution may be approximated by a non-interacting 
resonance gas obeying Boltzmann statistics. The pressure is then given by, 

^) ^/B(r)(cosh(^B/r)-i) , (6) 

/ res. gas 

where fs (T) contains all the information on the resonance mass spectrum. Within this 
approximation the ratio of Taylor expansion coefficients is given by (c„+2/c„)rcs. gas = 
9/((n + 2)(7i + 1)) and thus is indeed temperature independent. Moreover, the ratios 
are not affected by details of the hadronic mass spectrum and, in particular, are not 
affected by the quark mass dependence of the hadron spectrum. A straightforward 
comparison with lattice calculations which generally are performed at unphysically 
large quark mass values thus is justified. Similarly, the ratio Ap/T'^Xg = (1 ~ 
1/ cosh(3^q/T))/9 is independent of the resonance spectrum /s(T) and does not vary 
with temperature for fixed fj^s/T. These estimates are shown as straight lines in 
Fig. ^ for temperatures below To. Of course, the simple form of the /is-dependence 
in the resonance gas, Eq. El does not lead to any critical behavior in its own. At 
all temperatures the quark number susceptibility will rise with increasing fiq/T and 
for fixed Hq/T it rapidly rises with temperature. Moreover, as in the resonance gas 
model the dependence of e.g. the pressure on fiq/T only appears through the factor 
cosh(3/ig/T) the radius of convergence of a Taylor expansion, of course, is infinite, i.e. 
lim {cn+2/ {cn)res. gas = 0. Any unambiguous evidence for the existence of a chiral 

n — >oo 

critical point in the QCD phase diagram thus should also show up as clear deviation 
from resonance gas behavior. 

A direct comparison of resonance gas model calculations for thermodynamic 
observables like the pressure or susceptibility with lattice calculations is somewhat less 
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Figure 5. Comparison of the resonance gas model prediction for the temperature 
dependence of the coefficient C2(T) with lattice results (left) and the resulting 
dependence of Xq/T'^ on temperature (right) for various values of the chemical 
potential. Shown are results from Ref. '"25; where the resonance gas is compared 
to lattice calculations up to 0{^iq) only. The different curves shown in the left 
hand figure correspond to A = 0.9, 1.0, 1.1 and 1.2 (top to bottom). 



stringent as it is necessary to adjust the resonance spectrum to the still unphysical 
conditions that, at present, could be realized on the lattice. The Taylor series 
discussed in section 3 is based on lattice calculations performed with quark masses that 
correspond to a pseudo-scalar (pion) mass of about 770 MeV. In this case the mass 
of the lightest baryonic state is about twice as large as the experimental value for the 
nucleon mass. This quark mass dependence of the baryon spectrum, mjy (mpg), has 
been modeled using the ansatz mH{mpg)/mH = 1 + A{mpg /niH)'^ |25| . where 
is the experimental value of baryon resonances and mpg is the pseudo-scalar meson 
mass used in the lattice calculations. A comparison with the first expansion coefficient 
which only depends on the baryonic part of the spectrum*, C2(T) = 9/B(r)/2, shows 
that the temperature dependence for T<Tq is well reproduced with A ~ 1. With this 
also the temperature dependence of the quark number susceptibility is well described 
by the resonance gas in the low temperature phase (see Fig.EJ. This comparison can 
be extended to the meson sector contributing to p/T^ at /is = or to other quantum 
number channels controlled by other combinations of the quark chemical potentials, 
for instance the iso- vector channel with /// = — Also in these cases lattice results 
for thermodynamic observables at low temperature are well described by a hadronic 
resonance gas model |15) . 



5. Conclusions 



We have shown that the thermodynamics of QCD at small values of the baryon 
chemical potential can well be studied in terms of a Taylor expansion of the 
logarithm of the QCD partition function. In fact, a 6*'* order expansion seems to 

* The leading term, co(T), receives contributions from the mesonic as well as the baryonic part of 
the spectrum. 
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be quite sufficient for a reliable determination of the change in pressure as well as a 
determination of the temperature dependence of the net baryon number density for 
A'-B^3T. The good agreement between lattice results in the low temperature hadronic 
phase and resonance gas models is on the one hand conceptually appealing as the 
resonance gas model has also been so successful in the description of experimental 
data on observed particle ratios in heavy ion collisions j2S|- On the other hand it 
shows that higher order calculations with high accuracy and/or calculations at smaller 
quark masses are needed in order to find deviations from resonance gas behavior and to 
deduce evidence for the existence of a chiral critical point in the QCD phase diagram. 
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